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We report ^^^'^^^Sb nuclear quadrupole resonance (NQR) measurements in the filled- 
skutterudite superconductor PrOs4Sbi2 in the temperature range of 0.05-30 K. The electric 
field gradients (EFG), Vzz and Vxx — Vyy, at the Sb site exhibit unusual temperature depen- 
dence below 30 K. To explain these features, we discuss the coupling between the Sb nuclear 
quadrupole moment and Pr 4/^-derived multipole moments. The observed T dependence 
of EFG is well explained by the CEF quasi-quartet consisted of Fi singlet and F4 triplet 
states in the cubic point group T\. These results, in turn, are indicative of the importance 
of the coupling between the ^^^'^^"^Sb quadrupole moments and the hexadecapole moment 
caused by the quasi-quartet. 
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Much interest is drawn to a rich variety of the unconventional nature of the superconduct- 
ing state in the Pr-based heavy-fermion superconductor PrOs4Sbi2 with Tc = 1.85 K.^ An 
intriguing feature of the fundamental physics of this compound is strong electron correlations 
attributed to the Pr 4/^ electrons, which is evidenced by both the large electronic specific 
heat coefficient 7 = 350 ~ 700 mJ/mol K^^'^ and enhanced cyclotron-effective masses.^ Since 
a large specific heat jump is observed at Tc, the heavy quasiparticles themselves are expected 
to take part in the formation of the cooper pair. 

Despite such a heavy-fermion character, the magnetic properties of this material have 
been well explained by a crystal electric field (CEF) model for single- ion Pr 4/2 states. The 
CEF level scheme in the Th point group has been established as follows:^"^ a nonmagnetic 
Fi singlet ground state (OK), a F^^^ triplet first excited state (8K), F^^^ triplet (135K) and 
F23 doublet (205K) higher energy excited states. Recent experiments have shown that this 
low-lying singlet-triplet scheme leads to a field induced antiferro-quadrupole ordered state at 

* E-mail address:tou@hiroshima-u. ac.jp 

^Present address: Faculty of Science, University of the Ryukyus, Okinawa 903-0213 



1/7 



J. Phys. Soc. Jpn. Letter 

high fields after the superconductivity disappears."^' Thus, the most foundamental questions 
for PrOs4Sbi2 are the properties of these two phases and their mutual relations. 

(2) 

It has been believed that the low-lying Ti singlet and T\ triplet states are the key for con- 
trolling the exotic ground state properties. From a group-theoretical perspective, Shiina has 
discussed the roles of multipolc degree of freedom arising from the singlet-triplet states which 
are regarded as a pseudo-quartet; this quartet causes the field induced ordered state. ^ From 
different viewpoints, Otsuki and co-workers have investigated the origin of the CEF splitting 
by taking account of p — / hybridizations;^ they also pointed out the low-lying singlet-triplet 
states are regarded as a pseudo-quartet leading to multipole orders. Furthermore, Matsumoto 
and Koga have argued that the pseudo-quartet state is responsible for the appearance of 
the exotic superconductivity.^'' In this context, it is essential to unravel the relationship be- 
tween the multipole degree of freedom caused by the pseudo-quartet and the unconventional 
superconductivity in PrOs4Sbi2. 

To probe the effect of multipole interactions of / electrons, nuclear magnetic resonance 
(NMR) and nuclear quadrupole resonance (NQR) are effective tools, since nuclear moments 
of ligand atoms around a rare-earth site can couple with multipole moments of / electrons 
with the same time-reversal. ^^"^^ In practice, such a coupling have been discussed in several / 
electron systems CeBg,^^ Np02,^^ and URu2Si2."'^'^ For PrOs4Sbi2, an importance of a coupling 
between Sb nuclear quadrupole moments and Pr 4/ multipole moments has already been 
pointed out in an Sb NQR study based on the fact that the temperature (T) dependences of the 
electric field gradient (EFG) and nuclear spin-lattice relaxation rate (1/Ti) are anomalous. 

In this paper, we report a detailed study of T dependence of the EFG at the Sb site in 
PrOs4Sbi2 using the Sb NQR. The Sb NQR measurements have been extended to all the 
transitions for ^^^Sb and ^^^Sb NQR spectra and to the T range of 0.05-30 K. We have found 
that the anomalous T dependence of EFG parameters can be explained by taking account 
of a singlet-triplet CEF state forming a quasi-quartet; the coupling between ^^^'^^^Sb nuclear 
quadrupole moments and the 4/^-derived hexadecapole moment caused by the quasi-quartet 
is important at low temperatures below 15 K.^ 

Single crystals of PrOs4Sbi2 were grown by an Sb-self flux method reported by elsewhere.^ 
Tc of the present sample is 1.83 K. The residual resistivity was obtained to be po < 7/if2cm, 
and the observation of the de Haas-van Alphen (dHvA) signal provides the mean free path of 
I ~ 2000 A,^ indicating the high quality of the sample. Some of the single crystals were crushed 
into grains/powder to improve the NQR detection sensitivity. In addition, two single crystals 
were also used for the present study; sample No.l has dimensions of L([110]) x pr([110]) x 
H{[mi]) 3.2 X 3.2 X 2.4 mm^ with a trigonal shape; sample No.2 has L([100]) x I^([010]) x 
if([001]) ~ 1 X 1 X 3 mm^ with a rectangular shape. ^^iSb (/ = 5/2) and ^'^^Sh (I = 7/2) NQR 
measurements were carried out on a crushed-grain sample and single crystalline samples using 
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a conventional pulsed spectrometer. The Sb NQR spectra were obtained by tracing the spin- 
echo intensity as a function of the frequency. The NQR measurements in the T range of 0.4-30 
K were carried out at Hiroshima University, while those in the T range of 0.05-1.5 K were 
done at Osaka University. As long as we use the crushed-grain sample, the NQR linewidth is 
not different from that of the single crystals; while it becomes approximately double when we 
use a fine powder. 

Figure 1 shows the -^^^'^^^Sb-NQR spectra for two isotopes in the paramagnetic state at 
4.2 K. The observed ^^^Sb NQR spectrum consists of two narrow lines around the resonance 
frequencies of ^"^^vi = 53.309 MHz and ^'^^V2 = 85.0930 MHz, while the ^^Ssb NQR spectrum 
consists of three lines around ^'^^vi = 42.161 MHz, ^'^^V2 = 49.910 MHz and ^'^"'u?, = 78.704 
MHz. The obtained values are consistent with the previous report. Here it should be noted 
that the large value of an asymmetry parameter r] ~ 0.45^^ leads to mixing of pure eigenstates 

for = differing by Am = 2. Here, m represents the magnetic quantum number with 
respect to the z axis of the principal system. That is, typical eigenstates are given by a linear 
combination of ^m-, = C'm^'m- For those reasons, each NQR peak frequency is 

not equivalent to the integral multiplication of the nuclear quadrupole frequency uq 

Figure 2 depicts the T dependence of the resonance frequencies ^^^t'j {i = 1,2,3). Both 
^'^^v'2{T) and ^^^z^3(T) increase with decreasing T, being consistent with the previous report. 
On the other hand, it is noteworthy that ^^^i/i(T) exhibits opposite T dependence. ^'^^Ui 
{i = 1,2) arc also dependent of T (not shown). These features indicate that rj as well as vq 
changes with T. 

In order to deduce more accurate NQR parameters, the resonance frequencies for ^^^'^^^Sb 
nuclei of PrOs4Sbi2 were calculated by exact diagonalization of the nuclear spin Hamiltonian 



maximum EFG, where the principal axes of the EFG are defined according to the previous 
paper.^^ The calculations yield ^^^vq = 44.177 MHz, ^^^uq = 26.817 MHz, and r] = 0.4586 at 
T = 4.2 K. To obtain the T dependence of rj and uq, we calculated the resonance frequencies 
for ^^^'^^^Sb at various temperatures using Eq.(l). The T dependence of ^^^i'q{T) and r]{T) is 
shown in Fig. 3. The remarkable feature in the NQR parameters is the strong T dependence 
of both ^^^i'q{T) and r]{T); it is the first case in our knowledge. 

The NQR frequency vqiT) in metals without / electrons shows normally a monotonic 
increase with decreasing T. In such a case, vqiT) is often discussed in terms of an empirical 
relation of yQ{T) = z/q(0)(1 — aT^/^) with a > O}^ Since this relation shows a monotonic 
increase with decreasing T from high T as high as a room temperature, it cannot explain 



matrix of ^^^'-'^^'^Sb as. 
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the abrupt increase of vq{T) below 10 K in PrOs4Sbi2. Alternatively, the T dependence of 
both i'q{T) and ?7(T) in PrOs4Sbi2 could be attributed to the coupling between the ^^^'^^"^Sb 
nuclear quadrupole moments and Pr 4/^-derived multipole moments. 

To discuss the coupling of the ^^^'^^^Sb nuclear quadrupole moments and Pr 4/^-derived 
multipole moments, we need a CEF Hamiltonian for the cubic point group in PrOs4Sbi2, 
which is given as Hcef = ^4(0° + 50|) + Afs{Ol - 210^) + A'q{01 - Of), where A^, A(^, 
and A'q are the CEF parameters and O™ the Stevens operators associated with the multipole 
moments. 

(2) 

In PrOs4Sbi2, the level splitting between the singlet Fi state and triplet F4 first excited 
states is estimated as A ~ 8K.*^ Since the excitation energies to higher levels are considered 

(2) 

as more than lOOK, the low T properties of PrOs4Sbi2 would be governed by the Fi-F^ 
pseudo-quartet state.®' ^ Such a quartet state has fifteen multipole moments consisting of 
three dipoles, five quadrupoles, three octupoles, and four hexadecapoles (See Table II in 
Ref.(8)). Those multipoles of / electrons can couple with -"^^-"^'^^'^Sb nuclear moments with the 
same time reversal. In the paramagnetic phase, only a hexadecapole that belongs to a 
Fi representation has a finite thermal average, so that it can contribute to a shift of the 
asymmetry parameter as well as the NQR frequency. In fact, it is not difficult to make up an 
invariant hyperfine coupling between and the ^^^'^^^Sb nuclear quadrupole moments on 
the basis of the Pr-Sb local bond symmetry.®' "^^""^^ 

Let us examine the effect of a hexadecapole in detail. The symmetry analysis allows 
us to write the EFG parameters with the following expressions, 

Fz^(r) = F^.(0) + a(lf°), 

V,,^yy{T) = V,,-yy{0) + b{H^), (2) 

making use of ^22(0), Vxx~yy{0), a, and b as fitting parameters. Here Vxx-yy is defined by 
Vxx-yy = Vxx — Vyy • Within the pseudo-quartet model, the thermal average of is given by 
a simple distribution between the ground singlet and the excited triplet as 

l+3e-AABT' (3) 
where fc^ is the Boltzmann factor and A the level splitting (A = 8K). In the following we 
shall compare this formula with the observed NQR results. 

Figure 4 represents the T dependence of Vzz and Vxx-yy, where we use the values of Sb 
quadrupole moments, ^^ig = -0.597 x lO'^® m^, and ^^Sg = -0.762 x lO'^® (m^).!® The 
solid lines arc calculated curves using Eqs.(2) and (3) below 15 K; the fitting parameters are 
obtained as a = 4.76x10^^, b = -5.26 xlO^^, V^ziO) = -2.042x10^2, Vxx-yy{0) = -9.34x102^. 

Clearly, the agreements between data and Vzz{T) and Vxx-yy (T) are very good below 10 
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K. Our data, therefore, provide evidence tliat the local electronic state at the Sb site at low 
temperatures is dominated by the singlet-triplet CEF states. Alternatively, the microscopic 
origin of the T dependence of the EFG parameters can be explained by taking account of the 
coupling of the ^^^'^^^Sb nuclear quadrupole moments with the hexadecapole moment under 
the same time reversal symmetry. 

The deviations above 15 K cannot be explained by the inclusion of the triplet F^^^ second 
excited state (135 K). Recently Goto and co-workers have detected a thermally activated 
rattling motion due to the off-center Pr-ion by means of ultra sonic measurements, where 
the rattling motion is not random, but has F23 symmetry.^^ Such an anisotropic motion 
might deform the 4/ electron charge distrbution, hence it would cause the decoupling between 
the Pr 4f hexadecapole moment and Sb nuclear quadrupole moment. For that reason, we 
assume the fitting parameters a and b in Eq.(2) are replaced as a{l — a/{l + w^r^)} and 
6{1 — /3/(l -|- w^r^)}, respectively, by introducing a decoupling term with the Debye type 
dispersion. Here a and /? are decoupling parameters, uj is the angular frequencies, r is the 
correlation time r = Toexp{E/kBT) with an attempt time tq and an activation energy E for 
a motion. The calculated curves (dashed lines) in the Fig. 4 reproduces well the experimental 
data; the fitting parameters for the rattling motion are obtained as a = 0.15, /? = 0.4, and 
£^ = 80 K, assuming tq = 8.8 x 10~^^ s^^ and uj/2-k = 42 MHz. Good agreements suggest 
that the rattling motion of Pr-ion is a possible origin for the deviations above 15 K. Since this 
thermally activated rattling motion quenches at low temperature below 15 K, it is evident 
that the low temperatiuc properties of this compound is governed by the multipole degree of 
freedom. Here, it is worth emphasizing that the present Sb NQR results give evidence for the 
singlet-triplet CEF level scheme reported previously the other CEF ground states cannot 
reproduce the present Sb NQR results. 

The present Sb NQR results also give evidence for the multipole degree of freedom caused 
by the CEF quasi-quartet state from microscopic viewpoints. The multipole degree of free- 
dom possibly leads to the multichannel Kondo effect that conduction electrons interact with 
/-electrons through orbital channels such as multipoles.^^ Such a multichannel Kondo ef- 
fect in the periodic lattice, i.e., periodic multichannel Kondo effect, might cause the exotic 
ground states controlled by magnetic fields in PrOs4Sbi2, i.e., the unconventional heavy- 
fermion superconductivity and the field induced order. In this context, multipole fluctuations, 
e.g., quadrupole fluctuations, arising from the multipole degree of freedom are a candidate 
for the medium of the Cooper pairing in PrOs4Sbi2. Actually, Kotegawa et. al. have pointed 
out that Pr 4/^-derived quadrupole fluctuations are a key to understand the appearance 
of the unconventional superconductivity in PrOs4Sbi2 from the unusual T dependence of the 
NQR relaxation rate in the paramagnetic phase. To prove the multipole fluctuations, further 
measurements are required. 
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In summary, we have performed ^^^'""^^^Sb NQR measurements in the filled-skutterudite 
superconductor PrOs4Sbi2 in the T range of 0.05-30 K. The EFG parameters, Vzz and Vxx — 
Vyy, at the Sb site exhibit unusual T dependence befow 30 K. To explain this feature, we 
have discussed the coupling between the Sb nuclear quadrupole moment and Pr 4/^-derived 
multipole moments; the T dependence of the EFG parameters can be explained by the coupling 
of the Sb nuclear quadrupole moment and the Pr 4/^-derived hexadecapole moment which 
has a finite thermal average in the paramagnetic state. In the present case, it has turned 

(2) 

out to be identical to the thermal average for the Fi-F^ CEF pseudo-quartet state. The 
evidence for the multipole degree of freedom provides a crucial clue to clarify the intimate 
relationship between the field induced antiferro quadrupole order and the unconventional 
superconductivity in PrOs4Sbi2. Further experimental and theoretical efforts are needed to 
unravel the microscopic origin of the multipole order and unconventional superconductivity 
in PrOs4Sbi2. 
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Fig. 1. i2isb and ^^Ssb-NQR spectra in PrOs4Sbi2 at T = 4.2 K. 
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Fig. 2. Temperature dependence of the ^^^Sb-NQR resonance peaks in PrOs4Sbi2. 
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Fig. 3. Temperature dependence of (a) ^^^Sb-NQR frequency ^^^fq and (b) asymmetry parameter 77 
in PrOs4Sbi2. 
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Fig. 4. Temperature dependence of the EFG parameters, (a) Vzz and (b) V^x — Vyy in PrOs4Sbi2. 
The sohd Unes are fits by Eqs. (3) and (4). The dashed curves represent fits with the contribution 
from the rattling due to Pr-ion off-center motion. (See text) 
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